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Recently, it was shown that the coupling of center-of-mass and relative motion in atomic systems 
leads to inelastic confinement-induced resonances (ICIRs) [Phys. Rev. Lett. 109, 073201 (2012)]. 

In the present work, the possible occurrence of ICIRs in quantum dots is investigated. Particu¬ 
larly, electron-hole and electron-electron two-body systems with long-range Coulomb interaction 
are considered using the material parameters of GaAs. ICIRs are identified for the electron-hole 
system verifying the universal nature of the ICIR and, additionally, resonances due to the coupling 
of center-of-mass and relative motion are found also for the electron-electron system. In analogy to 
the coherent molecule formation appearing at ICIR in atomic systems a significant change in the 
mean distance between electrons and holes at the resonance is observed. By using the redistribution 
of the particle densities at the resonance position in modern quantum-dot experiments, the ICIR 
can provide a new technique for the control of the electron distribution in quantum dots and for the 
generation of single photons on demand. 


I. INTRODUCTION 

Artificial atoms and quantum dots (QDs) have at¬ 
tracted a great deal of attention during recent years [T]. 
Their modifiable physical properties offer a variety of new 
technologies in electronics and optoelectronics, such as 
single-photon sources 0131 or single-electron transistors 
0|5]. In QDs, the charge carriers, electrons and holes, 
are strongly confined in each spatial direction by an ex¬ 
ternal potential leading to a discretization of the energy 
spectrum. The external confinement can be realized by 
modern fabrication techniques like lithography and epi¬ 
taxy or electrostatically. 

While in theoretical treatments a simplified external 
potential such as a box potential 0 0 or a harmonic 
potential iHlo] is often adopted, the exact shape of the 
confinement of a quantum dot is in general not known. 
Yet, the confinement is definitely finite and hence never 
purely harmonic. Such an anharmonicity leads to a cou¬ 
pling of the center-of-mass (c.m.) and relative (rel.) mo¬ 
tion. It was demonstrated recently for ultracold atoms 
mma that the c.m.-rel. coupling allows for a controlled 
transfer of two atoms into a bound molecular state. This 
phenomenon is denoted as inelastic confinement-induced 
resonance (ICIR) [13]. 

In the present work the appearance of resonances due 
to c.m.-rel. coupling and thus ICIRs for a two-body sys¬ 
tem of Coulomb-interacting particles is investigated. The 
occurrence of ICIR is demonstrated, thus further demon¬ 
strating their universality. It is also shown that the ab¬ 
sence of a well-defined last bound state as is the case 
for ultracold atoms leads to modifications reflecting the 
influence of the long-range Coulomb interaction. 

In fact, ICIRs occur both for attractive interactions 
(excitons) and for repulsive interactions (electron pairs). 
For electron pairs it is demonstrated that the com.-rel. 
coupling does not cause a significant change in the mean 
interparticle distance but nevertheless such a coupling 
can result in instabilities and modifications of the charge 
distribution. In the case of excitons, a variation of the 


QD geometry allows, however, for a substantial modifi¬ 
cation of the mean distance between the electron and the 
hole. It is demonstrated that an adiabatic transition from 
an unbound electron-hole pair to a more tightly bound 
state seems experimentally feasible in an exciton system. 
Since a reduced mean distance changes the recombina¬ 
tion rate [14] , ICIR are expected to deliver a novel kind of 
deterministic single-photon sources. Excitons in electro¬ 
static traps might be especially suitable for this purpose, 
since they can be manipulated in situ. 

The paper is organized in the following way. In Sec. |TI| 
the mechanism of ICIR is introduced briefly for ultracold 
atoms. In Sec. [TIT] the full two-body energy spectrum of 
a confined Coulomb system and the computational ap¬ 
proach are introduced. Sec. |IV| describes the model for 
the QD confinement. The results are presented in Sec. 
Y\ T his includes the discussion of the energy spectra in 
V A[ the analysis of the coupling strength in Sec. |VB[ 


and the consequences of a variation of the confinement 
in Sec. V C| Fina lly, the paper closes with a conclusion 
and outlook (Sec.jV^. 


II. INELASTIC CONFINEMENT-INDUCED 
RESONANCES 

Recently, ICIR were discovered in the context of 
trapped ultracold atoms. Large losses of the trapped 
atoms for a given interatomic interaction strength (tuned 
with the aid of magnetic Feshbach resonances) were ob¬ 
served in m and explained by the occurrence of an ICIR 
in m- In a dedicated experiment [12] the occurrence of 
the ICIR was then more directly confirmed. 

The interatomic interaction in the ultracold regime is 
well described by the Fermi-Huang 6 pseudopotential 

TT ! \ 47rfi.2a . . d 

Us{r) = -(5(r) —r (1) 
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where m is the atomic mass and a the 5-wave scatter¬ 
ing length describing the collision in the limit of vanish- 
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ing collision energy (temperature). Thus the interparti¬ 
cle interaction strength is fully characterized by the sin¬ 
gle parameter a. Experimentally, it is possible to vary 
this parameter to almost arbitrary values using mag¬ 
netic Feshbach resonances m- The spectrum of two 
ultracold atoms confined in a harmonic trapping poten¬ 
tial that interact via the S pseudopotential is analyti¬ 
cally solvable [m HH] and is shown in Fig. for the 
example of an isotropic harmonic potential and for a 
variation of the inverse scattering length d^o/a where 
dho = \/d/{iiUj) is the harmonic-oscillator length with 
the harmonic-oscillator frequency uj and the reduced 
mass /i = 17111712 /{mi + m 2 ) of the two particles with 
masses mi and m 2 respectively. The reduced mass is 
equal to m/2 for the here considered case of atoms with 
equal masses. 




Figure 1. (Color online) Eigenenergy spectra of two atoms 
interacting via a S pseudopotential and confined in a 3D 
isotropic harmonic potential for varying (inverse) interaction 
strength represented by the scattering length a. The upper 
part shows the rel. motion spectrum, while the lower part 
shows the total spectrum including the c.m. motion. A non¬ 
vanishing coupling W between c.m. and rel. motion as in¬ 
troduced by, e.g., an anharmonicity of the external potential 
makes the crossings between states become avoided as illus¬ 
trated in the inset. 


A characteristic feature of the rel. motion spectrum 
is the occurrence of a shallow molecular bound state for 
positive values of the s-wave scattering length a (the state 
bending down to negative infinity in the rel. motion spec¬ 
trum in Fig. [^. If no c.m.-rel. coupling is present, as is 
the case for a harmonic confinement, the energy spec¬ 
trum of the full two-body system is obtained by adding 
the energies of the rel. and c.m. motion, respectively. As 
a consequence, molecular bound states with c.m. exci¬ 
tation cross with states of unbound atom pairs denoted 
as trap states, i. e. states above l.bhw in the rel. motion 
plot. In case of a vanishing c.m.-rel. coupling {W = 0) 
the states cross diabatically as indicated by the black 
dashed lines in Fig. and is visible in the inset in Fig. 
for W = 0. 

A coupling W 7 ^ 0, e. g. induced by an anharmonic 
trap, leads to an avoided crossing and thus allows for an 
adiabatic transition (red solid line in Fig. of the trap 
state into a molecular state. This transfer into the bound 



Figure 2. (Color online) Sketch of an avoided energy crossing 
of a molecular bound state with c.m. excitation and a trap 
state in the c.m. ground state. Passing through the crossing 
adiabatically (on the red solid lines) allows for an transforma¬ 
tion of the bound state into a trap state and vice versa. The 
black dashed lines indicate the diabatic curves. 


state and thus a two-body recombination is only possi¬ 
ble, because the excess binding energy can be transferred 
into c.m. excitation energy due to the anharmonicity of 
the external confining potential. This redistribution of 
binding energy and kinetic energy is an inelastic process 
and thus these c.m.-rel. coupling resonances are denoted 
as inelastic confinement-induced resonances [13]. It was 
demonstrated in [ 12 ] that the adiabatic transfer into the 
molecular bound state and the resulting reduction of the 
interparticle distance by one order of magnitude can be 
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performed fully coherently and controlled by tuning the 
scattering length using a magnetic Feshbach resonance. 

An alternative way to reach the resonance is to vary 
dho, i- e. the geometry of the external confinement which 
is feasible in Coulomb systems, e. g., for excitons in elec¬ 
trostatic traps. In contrast to atomic systems, where 
in the rel. motion spectrum of the S pseudopotential 
only a single bound state exists, attractively interact¬ 
ing Coulomb systems like, e. g., excitons consisting of an 
electron and a hole, possess an infinite number of bound 
states that have the character of the ones of a hydrogen 
atom, i. e., they show a Rydberg series in free space ap¬ 
proaching the continuum threshold. Inserting an exciton 
in a trap potential breaks the Rydberg series and leads to 
a smooth transition of bound to trap states. Hence the 
clear distinction of bound and trap states as it is present 
for the case of ultracold atoms is not possible. For re¬ 
pulsively interacting Coulomb systems, e. g., an electron- 
electron pair, there are no bound states present at all. 
A discretization is only induced by the trap potential. 
Thus, the different structure of the spectra of QDs pro¬ 
vokes the natural question whether ICIR occur in a QD 
and if an analog to the molecule formation in the atomic 
system can be observed. 


III. COMPUTATIONAL APPROACH FOR THE 
ENERGY SPECTRUM OF THE COULOMB 
SYSTEM 


A system of two interacting particles with the absolute 
coordinates ri and r 2 is described by the Hamiltonian 


=Ti(ri) +T2(r2) + Vi(ri) + V'2(r2) 

+ C/(|ri-r 2 |), (2) 


where Ti, T 2 , Vi, and V 2 denote the kinetic and poten¬ 
tial energies of the particles, respectively. The latter two 
represent the confinement of the QD. U denotes the in¬ 
terparticle interaction. Since the Coulomb interaction 


^Coul(^) = ± 


1 

67 - 


47reo|ri - r2| 


(3) 


with er as the dielectric constant depends solely on the 
interparticle distance, the Hamiltonian is expressed in rel. 
and c.m. coordinates, r = ri — r 2 and R = |(ri r 2 ), 
respectively. 


H{r, R) = JYrei(r) + i^cm(R) + IF(r, R) (4) 

iJrel(r) = Trel(r) + Fre^r) + Ucoul{r) (5) 

ifem(R) = Tc„,(R) + Vm(R) • ( 6 ) 


Here, the potentials Hrei(i*) and Hcm(R) represent the 
separable parts of the external potential, whereas the 
coupling term IF(r, R) includes all non-separable parts. 
Thus the Hamiltonians H^ei and Hem correspond effec¬ 
tively to single-particle systems and will occasionally be 
denoted in the following as single-particle Hamiltonians. 


For particles with identical mass confined to a har¬ 
monic potential the coupling term W (r, R) vanishes and 
the Hamiltonian is composed of the rel. and c.m. mo¬ 
tion of the two decoupled single-particle Hamiltonians, 
respectively. While the solution of the c.m. motion is 
well-known for a harmonic confinement, due to the long- 
range behavior of the Coulomb interaction the rel. mo¬ 
tion possesses only exact analytic solutions for certain 
energy levels using particular values of confinement and 
interaction strength [iQiiiaiio]. 

In this work, the eigenenergies and wave functions of 
the stationary Schrodinger equation with the Hamilto¬ 
nian Q, 


H\^i)=Ei\^i), (7) 

are calculated using an exact diagonalization approach 
[m [22]. Herein, the solutions of c.m. and rel. motion, 
i. e. the decoupled parts of , with their respective wave 
functions '^(R) and 0(r) are calculated separately. These 
wave functions are expanded in B splines for the radial 
part and spherical harmonics for the angular parts. The 
product states 4>,^(R, r) = (R) (r) form the basis 

of the solution of the full Hamiltonian. Thus, the full six¬ 
dimensional wave functions of the two-particle system 

«'i(R,r) = ^Ci,«$,(R,r) ( 8 ) 

Hi 

are constructed as superpositions of the In or¬ 

der treat computationally efficiently not only single 
wells (isolated quantum dots or ultracold atoms in a 
single-well potential), but also quantum-dot molecules 
or ultracold atoms in optical lattices, the basis func¬ 
tions are symmetry adapted to the eight irreducible 
representations of the orthorhombic point group D 2 h. 
(^Ag^ B\g^ B 2 g^ B^g^ B 2 U') B^u). This leads to a 

corresponding block structure of the Hamiltonian matrix. 

For locating ICIRs, a distinction of the different kinds 
of bound states for the Coulomb system is necessary. 
Similarly to the atomic system the mean interparticle 
distance r can be chosen as a measure of the binding 
strength for the attractively interacting electron-hole sys¬ 
tem. Thus, in order to characterize the Coulomb states, 
the mean radial distance 

00 

drrp(r) (9) 



is considered where 

p(r) = r2 y dVad^^r|4'(r,R)P (10) 

is the radial pair density. The dUa denotes the c.m. vol¬ 
ume element and dUr the angular volume element of the 
rel. motion. A state is regarded as a bound state, if f is of 
the same order of magnitude as the effective Bohr radius 
I (a^/'^o) where ttiq is the rest mass of an elec¬ 
tron and as the Bohr radius. If r is significantly larger 
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than a^, one can consider this state as weakly bound and 
its behavior is expected to be dominantly determined by 
the external trap potential. 

As introduction into the energy spectrum of Coulomb 
systems, two Coulomb-interacting particles of equal mass 
within the harmonic potential 

l/i(ri)= ^ Vjk]f (11) 

j=x,y,z 

are considered. With the aid of the parameters Vj and k 
the potential depth and the size of the QDs are adjusted 
to the desired QD confinement. A characteristic length 
scale of the potential is given by the harmonic-oscillator 
length dj = y^H/{jJioJj) defined for each spatial direction 
j = x^y^z with the harmonic-oscillator frequency uoj = 
kj along direction j. For an isotropic harmonic 
confinement for which V = Vj and k = kj and thus 
dho = dj and uj = ujj applies, the energy spectrum of the 
two Coulomb-interacting particles is shown in Fig. 

Only wave functions of the total-symmetric irreducible 
representation ag of the decoupled single-particle spec¬ 
tra, shown in Fig. |^a) and (b), are adopted when con¬ 
structing the energy spectrum of the full six-dimensional 
two-particle system in Fig. §0). InFig.lfd) the com¬ 
plete two-particle energy spectrum of Ag symmetry is 
constructed by using the single-particle wave functions of 
all symmetries. The orthorhombic D 2 h symmetry is only 
a subgroup of the proper symmetry group of the here con¬ 
sidered spherical-symmetrical problem that results from 
the isotropic harmonic trap and isotropic interparticle 
interaction. The orthorhombic group is chosen since in 
this work mainly anisotropic systems are considered and, 
as discussed earlier, the code allows for the treatment of 
generically orthorhombic problems like atoms in optical 
lattices or quantum-dot arrays. 

The parameter a = d/a^ can be positive or negative 
depending on whether the Coulomb interaction is attrac¬ 
tive or repulsive, representing an exciton system (nega¬ 
tive (t) or an electron pair (positive a ), respectively. The 
value of CF depends on included in . A change in 
can be considered as a screening effect. A variation of the 
interaction strength is a result of the modified properties 
of the environment, such as a modified charge density 
in the QD. However, a complete screening caused by , 
which suppresses the long-range behavior of the Coulomb 
potential, is usually not obtained in semiconductors [23] . 

The c.m. spectrum in Fig.[^(a) reveals the well-known 
simple harmonic-oscillator solution 

^ho = (2n + / + 2^diu , (12) 

where the quantum numbers n and I are restricted to 
0 < n, / (/ with only even numbers in ag symmetry). The 
energy levels are not influenced by a change of the inter¬ 
action strength. The relative-motion spectrum shown in 
Fig. I^b) varies with the interaction strength. At cr = 0, 
the Coulomb interaction is zero and the harmonic oscil¬ 
lator energies are revealed. For the even ag symmetry. 
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Figure 3. (Color online) Eigenenergy spectra for an exci¬ 
ton (cr < 0) or an electron pair (cr > 0) with me,h = 
0.067 mo confined in a harmonic potential with the dimen¬ 
sions {dx,dy,dz) = (31.6, 31.6, 31.6) nm . The parameter cr is 
modified by a change in Cr contained in a^. In (a) and (b) the 
ag energy spectra of the rel. and c.m. motion, respectively, are 
shown. The Ag energy spectrum of the full Hamiltonian 0 
for VF = 0 is shown when adopting only single-particle wave 
functions of ag symmetry (c) or single-particle wave functions 
of all symmetries of the D 2 h group (d). 


the ground state with (n,/) = (0,0) appears. The next 
energy level at cr = 0 comprises the states with quan¬ 
tum number (n,/) equal to (1,0) and (0,2), and thus, 
it is fourfold degenerate. The degeneracy comes from 
the quantum number m of the spherical harmonics with 
but only even values occur due to the con¬ 
sidered ag symmetry. A ninefold degenerate state with 
(n,/) corresponding to (2,0),(1,2) and (0,4) follows, etc. 
For positive values of cr the repulsive Coulomb interac¬ 
tion pushes the energy levels to higher values until they 
pass the continuum limit. For the attractive Coulomb 
interaction, i. e. cr < 0, the states bend down in energy. 
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Thereby the states with I = 0 react more sensitively than 
higher I > 0 states because of the presence of the cen¬ 
trifugal barrier for / >0. This leads to the occurrence of 
unavoided (real) crossings of states with different values 
of /. In the limit in which the trap depth approaches 
zero, the attractive Coulomb interaction reproduces the 
Rydberg series again. 

The Ag energy spectrum of the full Hamiltonian ([^ 
shown in Fig. I^c) is only composed of ag rel. and c.m. 
states. Since the c.m.-rel. coupling vanishes for this sys¬ 
tem, the spectrum is obtained by adding the two single¬ 
particle spectra of the c.m. and rel. motion. If all sym¬ 
metries of D 2 h from the rel. and c.m. motion are used for 
the construction of the full energy spectrum, as shown 
in Fig.j^d), some additional states appear that are (ap¬ 
proximately) built from single-particles states of the same 
symmetry, e. g., big (g) big = Ag [22]. The lowest energy 
levels still belong to the rel. and c.m. states in ag sym¬ 
metry as the comparison with the spectrum in Fig. [^c) 
shows. 

The Ag symmetry is chosen, since it contains the 
ground state. If the coupling term W (r, R) does not van¬ 
ish, all symmetries should be included in the full energy- 
spectrum calculations since they contribute to the com¬ 
position of the full wave function ^^(R, r). 

IV. C.M.-REL. COUPLING 

In order to simulate a more realistic single-well con¬ 
finement of a QD, the sextic potential [21] 

nr)= E + (13) 

j=x,y,z ^ ^ 

is chosen. This single-well potential introduces the re¬ 
quired anharmonicity for c.m.-rel. coupling. In fact, the 
anharmonicity of this sextic potential is relatively weak 
and thus the results discussed in this work are expected 
to be a rather conservative estimate for the strength of 
the coupling between c.m. and rel. motion as occurring in 
many real quantum dots, especially those that are more 
appropriately modeled by a square-well potential. Since 
in the numerical approach the sextic potential is obtained 
by an expansion of a sin^ potential up to the sixth or¬ 
der, a parameter for the strength of the anharmonicity 
is the ratio Vj/{bwj). The deeper the potential, i. e. the 
larger Vjl{hujj) :$> 1, the smaller the c.m.-rel. coupling 
(for fixed masses of the particles). 

Due to the anharmonicity of the confinement, the cou¬ 
pling term IF(r, R) in Q does not vanish and is com¬ 
posed of non-separable parts of the form with 

n^m e N\{0}. If the two particles have equal effective 
masses, such as the electron pair, the non-separable part 
consists of a polynomial with even values of n and m, i. e. 

r‘jRj and VjR^j. In the case of an exciton, where 
usually the electron and the hole have different effective 
masses, odd values of n and m also appear in IF(r, R), 
i.e. VjRj, VjR^, 


The matrix element 

Wa,^ = («'(“)(r,R)|VF(r,R)|«'(«(r,R)) (14) 

defines the coupling strength between the states and 
In order to examine the states involved in a cou¬ 
pling at the resonance position, their wave function den¬ 
sities are explored. 


V. RESULTS 

In order to obtain ICIRs comparable with those in 
atomic systems [TT] [12], a strong anisotropic QD con¬ 
finement is favored. In this work, inspired by the work 
on atomic systems, a cigar-shaped potential is chosen. In 
this case the longitudinal and transversal excitations ex¬ 
perience different trap strengths and thus have distinct 
extensions in the spatial directions providing a signifi¬ 
cant difference in the mean distance for a transition be¬ 
tween different excitations of them. Since a very large 
anisotropy is numerically challenging for the adopted 
computational approach, the longitudinal size, i. e., the 
size in x direction, is only chosen up to V^-times larger 
than the transversal sizes. Such cigar-shaped potentials 
are encountered in quantum dashes [24], nanorods [25] . 
and in the quasi ID regime also known as quantum wires 

m- 

All energy spectra are calculated for the Ag symmetry 
and taking into account single-particle wave functions of 
all symmetries. In the case of excitons the particles usu¬ 
ally possess different masses and are anyhow distinguish¬ 
able. In the case of an electron pair the particles are 
indistinguishable Fermions. Therefore, states with Ag 
symmetry in the spatial part correspond to spin-singlet 
states in the case of electron pairs. 


A. Variation of the interaction strength 

In most cases, the effective mass of an electron is much 
smaller than the effective mass of the hole. However, first 
an electron and a hole in a sextic potential are consid¬ 
ered where the hole mass is assumed to be equal to the 
effective electron mass, rrih = rrie = 0.067mo- Here, the 
c.m.-rel. coupling effects can only be caused by the anhar¬ 
monicity of the confinement, and thus, a clear distinction 
from other effects is possible. 

The fully coupled energy spectrum with variation of 
the Coulomb interaction strength is shown in Fig. Qa). 
An avoided crossing can be identified in the framed box 
and is shown magnified in Fig. [^b). This avoided cross¬ 
ing, here labeled with'd, with a full width at half maxi¬ 
mum (FWHM) « 0.018 has the position « —1.78 

with ^ 17.5 nm. 

The cuts through the wave-function density of the cou¬ 
pled states shown in Fig. |^c) and (d) demonstrate a 
clear redistribution of the c.m. and rel. excitations. For 
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Figure 4. (Color online) The Ag eigenenergy spectrum 
of the Hamiltonian 0 for an exciton with rrih = me = 
0.067 mo confined in a sextic potential with the dimensions 
(dx, dy, dz) = (100, 31.6, 31.6) nm is shown in (a), and the 
avoided crossing between coupled states, labeled with d, is 
magnified in (b). Cuts through the wave-function density 
IT(xe, Ze^Xh^Vh, Zh)\‘^ of the couplcd states along the x di¬ 
rection in absolute coordinates, Xh and Xe {ye,h = Ze^h = 0), 
in (c) and along the transversal x direction in the absolute 
coordinates, Zh and Xe {xe,h = ye,h = 0), in (d) confirm the 
c.m.-rel. coupling. 


a further illustration, a mapping of the states on famil¬ 


iar looking states of the c.m. and rel. motion from the 
Hamiltonian 0 can be performed. Looking at the cuts 
of the wave-function density from Fig.|^c), second-order 
c.m. excitations in the transversal direction can be de¬ 
duced from the three maxima of the density along the 
diagonal with the unidirectional z coordinates. Thus, the 
c.m. motion can approximately be described by the no¬ 
tation '0n(R) with n = {ux^ Tiy^ Uz) = (0, 2, 0) or (0,0, 2) 
inspired by the quantum numbers n of the harmonic os¬ 
cillator. Excitations in rel. motion can be revealed by 
looking at the diagonals with counterpropagating coor¬ 
dinates, either in x or in ^ direction, respectively. For 
the state in Fig. Qc), there is no rel. excitation, i. e., it 
is the ground state 0o(i*). This bound state with a mean 
distance r ~ 16 nm makes an adiabatic transition into 
a weaker bound state with rel. excitation, but no c.m. 
excitation, i. e. '0(o,o,o)(I^) which has r ^ 59 nm. 

This second state involved in the avoided crossing can 
be clearly identified from the cuts of the wave-function 
density in Fig. [^d) revealing two rel. excitations in the 
longitudinal direction but no c.m. excitations. Both cou¬ 
pled states posses an even symmetry in the rel. and c.m. 
motions, respectively. This result is consistent with the 
condition that only states with the same symmetry can 
couple as is required by the even powers of r and R in 
the coupling term W (r, R) for equal-mass particles. This 
type of avoided crossing d corresponds thus to the ICIR 
found for atomic systems. 


Further ICIRs between the transverse c.m. excited 
state and higher longitudinal rel. excited states occur. 
Here, a larger modification of the mean distance could 
be achieved. However, couplings between highly excited 
states are much weaker or even negligible due to the much 
smaller value of the integral in (14) that expresses the 
coupling strength. 

In the following, the particle-mass values are adjusted 
to realistic material properties of GaAs [27] . The system 
of an electron and a hole with the mass rrih = 10 rrie = 
0.67 mo is investigated. Its energy spectrum is shown in 
Fig-i in which two avoided crossings are marked. 


The green dashed-framed avoided crossing presented 
in detail in Fig. is labeled with Xq and has a FWHM « 
0.059 at position d^/a^ ^ —0.90, where « 34.5 nm. 
Here, the weakly bound state with a mean distance of 
r ~ 80 nm couples to a bound state with r ~ 24 nm. 
A similar significant change in mean distance is found 
for the orange-framed avoided crossing labeled with Aq 
and shown in detail in Fig. This crossing has a 
FWHM 0.006 and is located at d^/a^ —1.05 with 

~ 30 nm. Here, the same bound state as involved in 
Ag goes over into a weak bound state with r 65 nm. 
The interesting feature of Aq is the zero probability of 
the event that the two particles are located at the same 
place, i. e. the node at zero. Hence, the recombination is 
blocked in the weakly bound state. An adiabatic transi¬ 
tion into this kind of states could have significant influ¬ 
ence on photon emission. However, taking into account 
the coupling strength, which is discussed in more detail 
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Figure 5. (Color online) The Ag eigenenergy spectrum of 
the Hamiltonian @ for an exciton with mu — 10 me = 
0.67mo confined in a sextic potential with the dimensions 
{dx,dy,dz) = (100, 31.6, 31.6) nm. The two squares indicate 
avoided crossings, which are shown in detail in Fig. (green 
dashed Aq) and Fig. (orange Ao). 


in Sec. |VB[ the avoided crossing Aq is more difficult to 
be realized than the first avoided crossing Aq- 


In general, the coupled states at both avoided crossings 
of Fig. are very different from the states of the system 
with equal-mass particles. The coupling terms that now 
contain also odd powers of r and R do not only cause 
avoided crossings between states, but also lead to a gen¬ 
eral change of the shape of the wave functions. Thus, a 
mapping of the states to comparable product states of the 
c.m. and rel. motion appears not to be straightforward. 


Looking at unequal-mass particles in a harmonic con¬ 
finement, avoided crossings can appear due to the non¬ 
vanishing coupling term VjRj. The coupling term again 
strongly influences the shape of the states. Here, only 
avoided crossings similar to Aq can appear. 

In addition to the electron-hole system discussed 
above, the electron pair in a sextic potential is considered. 
The energy spectrum shown in Fig. [^a) reveals avoided 
crossings. Due to the repulsive Coulomb interaction and 
the confinement the states change only slightly in the in¬ 
teraction strength. Hence, only coupling between states 
with different longitudinal excitations is observed. One of 
these avoided crossings is labeled with ^ and magnified in 
Fig. [^b). Here, no significant change in the mean inter¬ 
particle distance occurs. The first transversally excited 
states cross only with very highly lying longitudinally 
excited states where due to the oscillatory behavior the 
coupling term becomes vanishing small. For electrons in 
a stronger confinement a crossing between transversally 
excited and low longitudinally excited states may, how¬ 
ever, be realized. Therefore, considering a QD array a 
transversal coupling to another QD due to the tunnel ef¬ 
fect is feasible. Consider the situation of a longitudinal 
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Figure 6. (Color online) The avoided crossing Aq framed by 
the big green dashed square in Fig. [^magnified in (a). Cuts 
through the wave function density |T^e, 2/e, Ze^Xh^yn^ Zh)\‘^ of 
the coupled states are shown along the absolute x coordinates, 
Xh and Xe {ye,h = Ze,h = 0), m (b) and along the absolute z- 
coordinates, Zh and Ze {xe,h = ye,h = 0), in (c). The c.m.-rel. 
resonance causes a significant change in distance between the 
electron and the hole. 


state that is well localized in a single QD and thus tun¬ 
neling to a transversally coupled QD is suppressed, since 
there is practically no electron density close to the tun¬ 
nel barrier. If this state is adiabatically transferred to 
the transversally excited state, this state has a non-zero 
tunnel probability to the transversally coupled neighbor 
QD. Evidently, an in-situ transfer form one state to the 
other may then be used as a switch that induces a charge 
migration from one QD to another one. 


B. Analysis of the coupling strengths 

To discuss and compare the time-dependent behavior 
at the avoided crossings, the Landau-Zener theory [28L 
|3Q] is used. The probability of an adiabatic transition at 
the avoided crossing can be estimated depending on the 
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Figure 7. (Color online) The avoided crossing Ao orange¬ 
framed in Fig. magnified in (a). Cuts through the 
wave-function density |T(a:e, 2/e, Ze^Xh^Vh^ Zh)\^ of the coupled 
states are shown along the absolute x coordinates, Xh and Xe 
{ye,h — Ze,h = 0), in (b) and along the absolute z coordinates, 
Zh and Ze {xe,h = Ve.h = 0), in (c). The c.m.-rel. resonance 
causes a transition of a state with zero probability for de,h = 0 
and a bound state. 


Landau-Zener velocity, i. e., the change in parameter a 
per time da/dt = 'T’lz- After calculating the linear func¬ 
tions for the involved diabatic states and the coupling 
between them, the upper bound of vlz can be estimated 
at which the adiabatic transition is more likely than the 
diabatic transition. 

Applying the Landau-Zener theory to the avoided 
crossing of the electron-hole system with equal-mass 
particles in Fig. the velocity vi,z,^ ~ 0.9 GHz is found 
as upper bound Tor an adiabatic transition. Regarding 
the required change in a for passing this avoided cross¬ 
ing, dcr ~ 0.1, the transition has to be performed within 
^ 0.11ns. This transition time is in the window of 
the exciton lifetime which is of the order of r ~ 1 ns in 
GaAs j3TH33] . This result indicates the potential impact 
of an ICIRs in QDs. 

Considering the avoided crossings of the electron-hole 
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Figure 8. (Color online) The Ag eigenenergy spectrum of 
the Hamiltonian Q for an electron pair with me,i = me ,2 = 
0.067mo confined in a sextic potential with the dimensions 
(dx, dy, dz) = (100,31.6,31.6) nm is shown in (a), and the 
avoided crossing between coupled states, label with G is mag- 
nihed in (b). Even for a repulsive Coulomb interaction the 
c.m.-rel. coupling causes avoided crossings. 


system with rrih = lOrUe in Fig. the upper bounds 
'^lz,Ag ~ 2.9 GHz for avoided crossing Xq and 'r’LZ,Ao ^ 
34 MHz for avoided crossing Aq are obtained. Since 
'^lz,Ag is larger than different mass values of the 

particles enhance the coupling strength and an adiabatic 
transition is dominant even for shorter time scales. Here, 
the transition time dty^ > 0.01 ns with a required change 
dcr 0.15 is obtained. Contrary, observation of the 
avoided crossing of type Aq is expected to be difficult due 
to the two orders of magnitude smaller value of 'C’lz,Ao 
in comparison to vlz,Ag- only the transition time 

d^Ao ^ 0.9 ns with the change daxo ~ 0-03. 

For the avoided crossing ^ of the electron pair from 
Fig. the velocity vlz,^ = 0.15 GHz is obtained. This 
velocity is considerably higher for transitions between 
two longitudinally excited states of the electron pair. 
Here, a change dcrj^^ 0.02 with the transition time 
dt^ > 0.13 ns has to be achieved. 

In general, the observation of an ICIR in experiments 
strongly depends on the good resolution of the param¬ 
eter cr. For this purpose, the uncertainty of the dielec- 
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trie constant [34], of the effective mass, and of the trap 
confinement must be kept low. As mentioned earlier, 
it should on the other hand also be reminded that the 
sextic potential chosen in the present simulations is only 
weakly anharmonic. Trap potentials with higher degree 
of anharmonicity like square-well potentials are expected 
to provide stronger couplings and thus broader avoided 
crossings. 


C. Variation of the confinement 


The variation of the dielectric constant in a QD to 
change the Coulomb interaction strength between the 
particles is experimentally a complicated task. A more 
suitable parameter for modification is the confinement of 
the QD. In electrostatic QDs, a modification of the con¬ 
finement is possible by a variation of the applied voltage 
[35] . The intrinsic problem arising in connection with 
electrostatic traps is the exciton dissociation by applying 
an electric field. However, various studies report elec¬ 
trostatic traps for indirect excitons in coupled QDs [36]- 
[39] . Also for QDs fabricated by chemical and growth 
processes the size can be varied in a controlled way [40] . 
Hence, it is appropriate to investigate a confinement vari¬ 
ation for the electron-hole system. 

In Fig. the confinement is varied from an isotropic 
to a cigar-shaped potential with x being the longitudi¬ 
nal direction. The parameter 77 = uj_\_Iux indicates the 
anisotropy. The energy levels bend down towards lower 
values the more anisotropic the confinement becomes. An 
avoided crossing labeled with k can clearly be located at 
77 « 2.2 and has a FWHM « 0.120. Here, a bound 
state with a small anisotropy of 77 = 1.9 and a mean dis¬ 
tance r « 18 nm couples to a weaker bound state with 
f ~ 45 nm at 77 ~ 2.4 nm. The densities (Fig. [^c) and 
(d)) of the avoided crossing k reveal the redistribution of 
c.m. and rel. excitations, the typical characteristic of an 
ICIR. The change in mean distance is not as significant 
as in the case of the variation of the Coulomb interaction 
in Fig. since the anisotropy at the ICIR is smaller. 
However, the position of the resonance may be tuned by 
varying the confinement strength in the transversal direc¬ 
tions and thus leading to other changes in mean distance. 
Moreover, for a system of particles with different masses 
a resonance with zero probability for de^h = 0 similar to 
the ICIR Ao can occur. 

Using the Landau-Zener theory, the upper bound 
~ 0.8 GHz for avoided crossing hi is found. This 
Landau-Zener velocity is similar to vi.z^'d obtained from 
the variation of the Coulomb interaction. With the 
change dr]^ « 0.5 in the confinement variation, the tran¬ 
sition time dt^ > 0.63 ns is found which is still smaller 
than the usual exciton lifetimes [3TII33] . This feature 
makes ICIRs induced by a confinement variation accessi¬ 
ble to an observation in experiments. 
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Figure 9. (Color online) The Ag eigenenergy spectrum of the 
Hamiltonian for an exciton with mu — rUe = 0.067 ttiq 
confined in a sextic potential with the transversal dimensions 
dy^z = 31.6 nm and dielectric constant Cr = 13.15 giving a 
d^jay ~ —1.52 is shown (a), and an avoided crossing be¬ 
tween coupled states, labeled with /^, is magnified in (b). Cuts 
through the density |T(xe, 7/e, ^e, 7//i, of the coupled 

states along the x direction in absolute coordinates, Xh and Xe 
{ye,h = = 0), in (c) and along the transversal z direction 

in the absolute coordinates, Zh and Ze {xe,h = ye,h = 0), in 
(d) are shown. 
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VI. CONCLUSION AND OUTLOOK 


It is demonstrated that inelastic confinement-induced 
resonances, first described for trapped ultracold atoms 
iiiiiia, occur in the electron-hole Coulomb systems 
confined in a QD. Furthermore, the c.m.-rel. coupling 
leads also to other types of avoided crossings as, e. g., 
in electron-electron systems. In order to investigate a 
realistic QD, its confinement is approximated by an an- 
harmonic and anisotropic potential. For the electron-hole 
system, the ICIR causes a significant change in the in¬ 
terparticle distance, and thus can lead to a significant 
influence on the recombination rate. The change in dis¬ 
tance grows with the anisotropy of the confinement, since 
couplings between weakly bound longitudinally excited 
states and bound transversally excited states become 
possible. Thereby, the anharmonicity of the confine¬ 
ment induces the c.m.-rel. coupling. Hence, the coupling 
strength can be controlled by the confinement geome¬ 
try. Furthermore, looking at realistic exciton lifetimes a 
transition between states at an ICIR appears to be in an 
experimentally feasible time window. 

An interesting application of the ICIR in Coulomb sys¬ 
tems studied in the present work could be a novel kind of 
single-photon sources on demand based on excitons. In 
the case of excitons in self-organized QD the in situ vari¬ 
ation of the confinement, i. e. the change of the trap po¬ 
tential on a timescale shorter than the exciton lifetime, is 


usually impossible. However, the realization of excitons 
in electrostatic traps [37H39] or laser-induced traps m 
overcomes this limitation and thus allows for real-time 
reduction of the relative mean distance between electron 
and hole at an ICIR. The enhanced recombination prob¬ 
ability leads then to a controlled photon emission. How¬ 
ever, in order to stabilize the exciton against dissociation 
an additional electric held perpendicular to the plane of 
electrostatic confinement is usually adopted. The influ¬ 
ence of this alignment of the exciton for the ICIR will be 
subject of a future study. In the case of more than one ex¬ 
citon, the alignment leads to a dipole-dipole interaction 
between the excitons [39]. Noteworthy, the occurrence 
of ICIR in cold dipolar gases has recently theoretically 
been demonstrated in |42|. For electrons in QD arrays 
the ICIR may on the other hand be used as a possible 
switch that allows to turn-on or -off the tunnel current 
to a neighbor QD. 

Moreover, this work confirms the universal nature of 
ICIR as they are now proven to appear in the two very 
contrary systems of quantum dots and ultracold atoms. 
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